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1.  INTRODUCTION 


Wrap-around  fins  have  been  used  primarily  for  their  advantages  in  packaging  tube- 
launched  projectiles.  The  wrap-around  fin  conforms  to  the  cylindrical  shape  of  the  projectile 
while  in  the  launch  tube,  allowing  more  efficient  use  of  space.  Thus,  greater  numbers  of 
wrap-around  fin  projectiles  can  be  stored  in  the  same  space  as  fixed  fin  projectiles  designed 
to  deliver  the  same  payload  (Dahlke  1975;  Winchenbach  19S6). 

The  cylindrical  shape  of  the  wrap-around  fin  is  advantageous  for  packaging;  but  it  can 
also  be  compromising  to  the  dynamic  stability  of  the  projectile.  In  many  configurations  that 
wrap-around  fins  have  been  employed,  it  has  been  noted  that  the  roll  moment  coefficient 
may  change  in  magnitude  and  sign  as  the  Mach  number  varies  (Dahlke  1975;  Winchenbach 
19S6,Mermagen  1981).  During  the  course  of  flight  of  a  wrap-around  fin  projectile,  it  is  possi¬ 
ble  for  its  spin  rate  to  increase  or  decrease  more  than  once.  In  addition,  the  direction  of  spin 
may  change.  This  type  of  behavior  can  produce  poor  flight  dynamics.  In  order  to  design  dy¬ 
namically  stable  projectiles  employing  wrap-around  fins,  it  is  necessary  to  have  the  ability  to 
predict  the  roll  moment  coefficient  at  all  flight  conditions  for  the  full  trajectory.  Design  code 
methodology  is  inadequate  for  this  problem  (Dahlke  1990).  Results  from  computational  fluid 
dynamics  (CFD)  calculations  have  shown  promise.  Normal  force  coefficients  calculated  us¬ 
ing  inviscid  CFD  computations  have  compared  favorably  with  experimental  data.  However, 
roll  moment  coefficients  calculated  using  inviscid  CFD  computations  have  not  shown  good 
agreement  with  experimental  data  (Wardlaw  1987).  This  paper  presents  the  initial  results 
for  establishing  the  cai)abilitv  of  predicting  the  roll  moment  coefficient  for  a  projectile  with 
wrap-around  fins  llirough  viscous  computations  with  a  3-D  full  Navier-Siokcs  code. 

2.  ABOUT  THE  TEST  CASE 

The  experimental  data,  used  for  comparison  were  obtained  from  references  1,  6  and  7. 
The  reports  document  a  comprehensive  cfloit  to  experimentally  determine  how  changes  in 
geometry  affect  the  aerodynamic  forces  generated  by  wrap-around  fuis.  .-X  standard  wrap¬ 
around  fin  projectile  determined  by  the  Technical  Cooperation  Program  (  TTCP),  as  seen 
in  Figure  1,  was  used  as  the  basic  configuration.  A.  number  of  geometric  variations  to 
the  basic  configuration  wore  made.  The  aerodynamic  forces  of  each  configuration  were 
measured  and  documented.  I'hc  configuration  for  the  wrap-around  fin  projectile  modeled 
in  the  computation  was  derived  fr^^rn  the  standard  set  by  the  TTCP  (Dahlke  1975,  Dahlke 
1976).  Although  the  body  retains  llie  dimensions  of  the  standard  TTCP  configuration,  the 
fins  differ  slightly.  The  standard  TTCP  configuration  had  fins  with  symmetric  icnding  and 
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trailing  edge  bevels.  The  computational  model  had  blunt  leading  and  trailing  edges.  There  is 
a  dilTcrcnceof  45  degrees  between  the  root  and  tip  chord  in  the  standard  TTCP  configuration 
while  the  root  and  tip  chord  are  parallel  in  the  computational  model.  See  Figures  2  and  4 
for  a  visual  reference.  Another  difTerence  to  note  is  that  the  standard  TTCP  configuration 
had  boundary  layer  trips  on  the  body  and  the  fin  leading  edges  (Dahlke  1975,  Dahike  1976). 
The  boundary  layer  was  laminar  for  the  body  and  fins  for  the  computational  model. 

Flowfield  solutions  were  computed  for  a  series  of  Mach  numbers  ranging  from  1.3  to  3.0 
at  sea-level  atmosidieric  conditions.  The  freestream  conditions  were;  density  =  1.19SA.i//;r?^; 
static  temperature  =  21.6°C;  and  static  pre.ssurc  =  101.3^'Pfi.  The  Reynolds  number,  based 
on  projectile  length,  for  sea-level  conditions  ranged  from  30  to  69  million,  Flowfield  solutions 
were  also  computed  using  wind  tunnel  free  stream  conditions.  The  Reynolds  number  under 
wind  tunnel  conditions  varied  from  17  million  at  Mach  1.3  to  69  million  at  Mach  3.0.  The 
wind  tunnel  atmospheric  conditions  were  obtained  from  reference  7.  The  experimental  data 
were  obtained  from  four  dilTcreiil  sites:  the  McDonnell  Douglas  .Aerophysics  wind  tunnel; 
the  .\ruold  Engineering  Development  Center  (.\EDC);  N,AS.'\  Langley  Research  Center;  and 
the  Jot  Propulsion  Laltoratory  (.IPL).  The  experiments  conducted  by  McDonnell  Douglas, 
ALDC  and  NASA  Langley  were  wind  tunnel  tests  while  the  JPL  data  were  obtained  from 
a  free  flight  test  (Dahlke  1975,  Dahlke  1976).  The  models  used  in  the  McDonnell  Douglas, 
AEDC  and  NAS.'\  Langley  experiments  followed  the  TTCP  specifications.  The  models 
used  to  obtain  the  experimental  data  at  JPL  were  scaled  versions  of  those  used  for  the 
wind  tunnel  experiments.  These  models  did  not  have  boundary  layer  trips  on  the  leading 
edges  of  the  fins  (Dahlke  1975.  Dahlke  1976).  The  experimentally  obtained  roll  moment 
coellicients  have  been  compared  to  values  computed  from  flowfield  solutions  calculated  by 
the  HRL  ZONAL  .Navier-Stokes  code.  Wind  tunnel  data  were  available  comparing  the  TTCP 
standard  configuration  with  a  coiiliguration  that  had  blunt  leading  and  trailing  edges.  See 
Figure  3  for  this  configuration.  The.sr  data  will  be  discussed  later. 

3.  CODE  AND  BOUNDARY  CONDITIONS 

.M!  flowfield  solutiiuis  were  corniiuted  with  the  URL  Zonal  code.  'J’he  URL  Zonal  code  is 
a  3-D  full  .N'avier-Si  Tes  code  which  can  be  a|/plied  tcj  zonal  grid  to[)ologie.s,  J'lie  governing 
ecjualion.s  are  niiniei  icallv  integrated  by  an  exirlicil  tinic-marching  method  (Patel.  Sturek 
and  Smith  19S9,  I’.itel  and  Edge  1991). 

The  flowfield  solutions  were  viscous  flow  calculaticuis.  The  computationally  modeled 
sill  face  was  oiie-foui  tli  of  t  h.;  cnliie  project  ile  .surlacc.  The*  symmetry  of  the  projectile  allowed 
the  iinnleliiig  of  one  fin  and  apiiroxiinatelv  forty-five  degrees  of  body  surface  on  eitlier  side 


of  the  fin.  A  periodic  boundary  condition  was  written  to  take  advantage  of  this  symn 
This  helped  to  reduce  the  number  of  grid  cells  needed  for  the  computation,  but  ii  restricted 
the  angle  of  attack  to  zero  degrees.  The  calculation  of  the  roll  moment  coefficient  wr^s  the 
primary  goal  in  this  initial  effort.  For  this  reason,  the  flowfield  was  computed  for  the  body 
and  fins  only,  since  the  influence  of  the  recirculating  flow  at  the  base  on  the  roll  moment 
coefficient  is  expected  to  be  small.  Future  computations  which  include  flow  field  solutions 
at  the  projectile  base  will  hopefully  verify  the  accuracy  of  this  cissumption.  A  zero  gradient 
boundary  condition  was  used  at  the  trailing  edges  of  the  fins  for  the  downstream  boundary 
condition.  A  i*on-reflecting  boundary  condition  was  applied  to  the  outermost  grid  plane 
from  the  body  surface.  The  non-reflecting  boundary  condition  allowed  the  outermost  grid 
plane  to  be  placed  relatively  close  to  the  body,  which  reduced  the  number  of  computational 
points  needed  for  a  solution.  In  order  to  obtain  a  viscous  solution,  the  fin  and  body  surfaces 
were  modeled  with  a  no-slip  boundary  condition. 

zonal  approach  was  used  to  obtain  a  solution  for  the  flowPeld.  In  the  zonal  approach, 
as  implemented  in  tlie  BRL  Zonal  code,  overlapped  zones  share  at  least  one  grid  cell  in  a 
given  direction  with  an  adjacent  zone.  The  shared  grid  cells  have  identical  coordinates  in 
both  of  the  overlapped  zones.  Since  the  two  overlapped  zones  have  the  same  coordinates  for 
the  shared  grid  cells,  zonal  coupling  requires  only  the  transfer  of  information  from  the  field 
of  one  zone  to  the  boundary  of  the  other  zone  and  vice-versa.  No  interpolation  is  required. 
A  multi-zone  solution  is  obtained  by  performing  the  integration  of  the  governing  equations 
in  all  zones  and  then  exchanging  information  between  overlapping  zones  before  advancing 
to  the  next  iteration  (Patel,  Sturek  and  Smith  1989;  Patel  and  Edge  1991). 

.All  compn'.  '.lions  were  performed  on  a  Cray-2  supercomputer.  Some  flowfield  solutions 
were  compu*  d  on  the  TACOM  Cray-2  while  others  were  computed  on  the  Ciay-2  at  BRL. 
As  configured  for  this  case,  the  BRL  Zonal  code  required  forty-five  million  words  of  memory. 

4.  GRID 

The  zonal  approach  facilitated  the  building  of  the  computational  grid.  The  unswept 
wrap-around  fins  with  blunt  leading  edges  would  have  been  difficult  to  model  with  a  wrap¬ 
around  grid.  The  use  of  a  wrap-around  grid  would  have  resulted  in  large  and  rapid  variations 
of  the  metric  terms,  which  could  have  severely  degraded  the  quality  of  the  solution  (Patel, 
Sturek  and  Smith  19S9;  Patel  and  Edge  1991).  The  zonal  approach  allowed  the  accurate 
modelling  of  the  wrap-around  fin  projectile’s  geometry  while  retaining  a  smooth  continuous 
computational  mesh.  .A  cut-away  view  of  the  computational  grid  can  be  seen  in  Figure  5. 
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The  grid  was  designed  for  viscous  compulation.  It  was  highK  clustered  near  the  body 
and  fin  surfaces.  The  total  nurtuer  of  points  used  for  the  computation  wa.s  951, SSS.  The 
dimensions  for  each  of  the  seven  zones  were  as  follows,  (20  x  80  x  32),  (20  x  SO  x  32),  (130 
X  80  X  32),  (130  X  80  X  32),  (20  x  80  x  18),  (96  x  SO  x  18),  (36  x  26  x  48).  .As  stated  earlier, 
these  seven  zones  were  configured  to  model  one-fourth  of  the  wrap-around  fin  projectile. 

An  algebraic  grid  generator,  developed  at  BRL,  was  used  to  build  the  computational 
mesh.  Once  the  interior  mesh  of  a  zone  is  completed,  grid  points  on  tlie  boundary  are 
added,  deleted  or  changed  to  create  overlaps  between  adjacent  zones.  In  order  for  the 
periodic  boundary  condition  to  function  pr  perly,  the  planes  in  which  data  were  exchanged 
were  spaced  such  that  they  were  exactly  ninety  degrees  apart  in  the  circumferential  direction. 

5.  RESULTS 

Qualitatively,  the  flowfield  solutions  computed  by  the  BRL  Zonal  code  show  a  number 
of  interesting  features  of  wrap-around  fin  aerodynamics.  Figures  C,  7  and  8  show  surface 
pressure  contours  on  the  concave  and  convex  .sides  of  the  wrap-around  fins  at  different  Mach 
numbers.  These  figures  show  the  iniluence  of  shocks  generated  at  the  fin  leading  edge  on 
the  fin  surface.  In  addition,  the  figures  also  show  the  development  of  a  high  pressure  region 
near  the  leading  edge  of  the  concave  side  of  the  fin  as  the  Mach  number  increases. 

The  differences  in  the  surface  pressure  contours  on  the  convex  and  concave  sides  of  the 
fin  are  quite  evident.  Figure  9  is  a  view  of  the  wrap-around  fin  from  its  concave  side.  Figure 
10  is  a  view  from  the  convex  side.  Figure  11  shows  a  view'  between  the  fins.  Near  the  fin 
root,  the  shock  generated  by  the  fin  is  clearly  asymmetric.  Figure  12  shows  a  flowfield  plane 
away  from  the  body  that  intersects  the  fin  at  approximately  fifty  percent  of  its  height.  It 
should  be  noted  that  the  plane  is  not  at  a  constant  radial  distance  from  the  body  surface. 
The  shock  on  the  concave  and  convex  sides  of  the  fin  appear  to  be  very  similar  to  each  other 
at  a  distance  away  from  the  body.  Near  the  fin  leading  edge  at  the  root  of  the  fin,  the  surface 
pressure  on  the  concave  and  convex  sides  of  the  fin  appear  to  be  very  different  from  each 
other. 

Figure  13  shows  the  flowfield  at  various  axial  stations  along  the  fin  length.  These  contours 
give  a  good  indication  of  the  overall  structure  of  shocks  generated  by  a  wrap-around  fin 
projectile.  The  contours  also  show  the  ''focusing”  of  the  shock  on  the  concave  side  of  the 
fin.  This  is  most  ajjparent  near  the  leading  edge  of  the  fin.  Generally  speaking,  the  leading 
edge  is  the  area  where  there  is  the  greatest  difference  in  pressure  between  the  concave  and 
convex  sides  of  the  fin.  This  is  in  agreement  with  the  findings  ineiitioned  in  Beference  1. 
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A  comparison  of  the  computed  roll  moment  coefRcient  and  the  experimentally  obtained 
roll  moment  coefficient  shows  good  agreement  with  overall  trends.  Figure  14  is  a  graph 
showing  the  roll  moment  coefficient  versus  Mach  number  for  experimental  data  obtained 
at  JPL  and  computations  made  from  the  BRL  Zonal  code  solutions.  Figure  14  shows  that 
the  roll  moment  coefficients  computed  from  BRL  Zonal  code  flowfield  solutions  follow  the 
general  trends  of  the  experimental  data.  At  lower  Mach  numbers,  the  roll  moment  coefficient 
is  positive.  At  higher  Mach  numbers,  the  roll  moment  coefficient  is  negative.  A  positive  roll 
moment  coefficient  indicates  a  roll  direction  towards  the  fin’s  center  of  curvature.  The 
predicted  and  JPL  measured  roll  moment  coefficients  indicate  one  cross  over  point,  a  Mach 
number  at  which  the  roll  moment  coefficient  is  zero,  near  Mach  1.7  . 

It  is  expected  that  the  geometrical  differences  between  the  computational  model  and  the 
TTCP  standard  configuration  would  produce  differences  in  their  roll  moment  coefficients. 
Additional  wind  tunnel  data  are  presented  to  show  how  the  fin  bluntncss  affects  the  roll 
moment  coefficient.  Figure  15  is  a  roll  moment  coefficient  versus  Mach  number  plot  for 
two  wrap-around  fin  configurations  obtained  under  similar  conditions.  One  configuration 
is  the  standard  TTCP  configuration  while  the  second  configuration  is  a  modified  TTCP 
configuration  which  has  blunt  fin  edges.  Figure  3  shows  the  fin  geometry  of  this  second 
configuration.  The  data  are  presented  because  the  fin  leading  edge  bluntness  is  the  primary 
geometrical  difference  between  the  TTCP  standard  configuration  and  the  computational 
model.  As  can  be  seen  in  Figure  15,  the  roll  moment  coefficient  versus  Mach  number  curve 
of  the  blunt  finned  configuration  seems  to  retain  the  same  overall  characteristics  of  the  curve 
for  the  TTCP  standard  configuration.  The  differences  in  the  two  curves  indicate  that  a 
blunt  leading  edge  as  opposed  to  a  forty-five  degree  leading  edge  wrap-around  fin  will  have 
a  slightly  higher  cross-over  point  and,  subsonically  at  least,  have  a  greater  roll  moment 
coefficient. 

A  comparison  of  the  different  roll  moment  coefficients  of  a  standard  TTCP  configuration 
obtained  from  various  facilities  shows  some  scatter  in  the  data,  but  the  overall  trends  stated 
earlier,  remain  intact.  Figure  16  is  a  plot  of  roll  moment  coefficient  versus  Mach  number  for 
experimental  data  obtained  using  the  TTCP  standard  configuration.  As  stated  earlier,  the 
McDonnell  Douglas,  AEDC  and  Langley  data  were  obtained  from  wind  tunnel  experiments 
with  the  standard  TTCP  configuration.  The  JPL  frec-flight  data  were  obtained  using  a 
scaled  down  version  of  the  TTCP  standard  configuration.  The  model  used  at  JPL  did  not 
have  boundary  layer  trips  on  the  fin  leading  edges  while  the  models  used  at  McDonnell 
Douglas,  AEDC  and  Langley  did  have  fin  leading  edge  boundary  layer  trips  (Dahlke  1975, 
Dahlke  1976).  In  the  Dahlke  1975  and  1976  reports,  the  data  in  Figure  16  was  presented  with 
a  discussion  of  the  effects  of  Reynolds  number  on  roll  moment  coefficient.  The  experimental 
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data  indicate  that  the  Reynolds  number  may  have  an  efTect  on  the  roll  moment  coefficient 
and  Dahlke  speculated  the  Reynolds  number  and  the  lack  of  boundary  layer  trips  may  be 
responsible  for  the  difference  in  the  JPL  data  iind  the  other  experimental  data.  The  JPL 
data  predicts  the  cross-over  point  to  be  at  approximately  Mach  1.7  while  the  other  wind 
tunnel  data  indicates  Mach  1.  With  the  exception  of  the  differing  cross-over  points,  the  JPL 
data  agrees  well  with  the  other  wind  tunnel  data.  The  computed  roll  moment  coefficients, 
shown  in  Figure  14,  indicate  that  there  is  little  difference  in  the  roll  moment  coefficient  under 
wind  tunnel,  (Reynolds  number  17  to  23  million)  and  sea-level,  (Reynolds  numbers  30  to  69 
million)  conditions. 

6.  CONCLUDING  REMARKS 

Flowfield  solutions  of  a  wrap-around  fin  projectile  were  calculated  with  the  BRL  3-D 
Zonal  Navier-Stokes  code.  The  roll  moment  coefficient  was  calculated  from  the  flowfield 
solutions  and  compared  to  the  roll  moment  coefficient  obtained  from  experiment  measure¬ 
ments  for  a  similar  wrap-around  fin  projectile  configu'‘ation.  The  computed  roll  moment 
coefficients  predicted  the  same  trends  seen  in  the  experimentally  obtained  roll  moment  co¬ 
efficients  in  magnitude  and  sign,  (direction).  Experimental  data  obtained  from  McDonnell 
Douglas,  AEDC  and  NASA  Langley  wind  tunnels  indicated  that  the  cross-over  point  was  at 
Mach  1  while  the  JPL  free  flight  data  indicated  that  the  cross-over  point  was  greater  than 
Mach  1.  It  is  encouraging  that  the  computed  roll  moment  coefficient  cross-over  point  was 
close  to  the  cross-over  point  of  the  JPL  free  flight  data.  The  boundary  layer  was  laminar  on 
the  fins  in  the  computational  model.  Also,  the  model  used  for  the  JPL  experiments  did  not 
have  a  boundary  layer  trip  on  the  fin  leading  edge,  while  the  models  used  in  the  wind  tunnel 
experiments  did.  Another  point  to  note  is  that  the  JPL  experiments  had  the  lowest  Reynolds 
number.  This  raises  the  possibility  that  the  wraji-around  fins  of  the  JPL  model  may  have 
had  a  laminar  boundary  layer  at  low  Mach  numbers  which  would  give  this  experiment  an¬ 
other  feature  in  common  with  the  computation.  The  experimental  data  indicate  that  the 
roll  moment  coefficient  varies  with  the  Reynolds  number.  However,  the  data  do  not  estab¬ 
lish  any  trends  which  would  allow  one  to  slate  conclusively  whether  there  is  a  relationship 
between  the  Reynolds  number  and  roll  moment  coefficient  or  what  that  relationship  may  be. 
The  results  from  the  CFD  computations  do  not  indicate  that  the  roll  moment  coefficient  is 
Reynolds  number  dependent. 

The  roll  moment  coefficient  computations  have  shown  good  agreement,  in  terms  of  the 
magnitude  and  sign,  with  values  of  the  roll  moment  coefficient  obtained  experimentally. 
The  computational  technique,  as  outlined  in  this  report,  has  shown  great  promise  in  being  a 


reliable  method  for  predicting  the  roll  moment  coefficient.  However,  a  one-to-one  comparison 
of  computation  and  experiment  is  difficult  to  make.  As  this  is  an  ongoing  effort,  further  work 
.will  be  done  to  fully  establish  this  capability.  Future  plans  include:  1)  applying  the  F3D 
code  to  wrap-around  fin  configurations  (Sahu  1988);  2)  computing  turbulent  viscous  flowfield 
solutions;  and  3)  computing  flowfield  solutions  for  small  angles  of  attack. 
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Figure  1.  Wrap-aroan<i  Fin  Model  Geometr 


Figure  2.  TTCP  Standard  Wrap-around  Fin 


Figure  3.  Blunt  Leading  Edge  Wrap-around  Fin 


Figure  4.  Computational  Model  Wrap-around  Fin 
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Figure  5.  Computational  Grid 
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Figure  9.  Nol•maIi^ed  Pressure  Contours  on  Concave  Side  of  Fin  at  Ma.ch=2 


Figure  10.  Normalized  Pressure  Contours  on  Convex  Side  of  F'in  at.  Mach— 2 
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Figure  11.  Normalized  Pressure  Contours  Between  F'in.s  at  Mach— 2 


Normalized  Pressure  Contours  Between  Fins  at  Mach— ‘2.5 
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Figure  13.  Normalized  Pressure  Contours  at  Varioiif;  Axial _S t ations  at  Mach  =2.5 


z 

UJ 

u 

LL 

LL 

LU 

o 

o 

K 

Z 

LU 

o 


o 

Cl 


0.04  p 
0.03 
0  02 
0.01  L 

i 

0.00  k 

I 

-0.01  ^ 

-0.02  p 

-0.03 

I 

I 

-0.04  - 


Rn  X  10”  Based  on  Length 

•  30-69  ZONAL  Code 
Xp  O  17-23  ZONAL  Code 


\  I 

V 

■'* 

'a 


2  5  JPL 


0.0  0.5 


1.0  1.5  2.0  2.5 

MACH  NUMBER 

Figure  14.  Roll  Moment  Coefficient  vs.  Mach  T'uiiiber  for 
Computational  and  JPL  Experimental  Data. 


3 . 5 


1- 

Z 

CJ 

LU 

U. 

LU 

O 

u 

H 

Z 

LU 

2 

o 

_J 

_l 

O 

cc 


Figure  15. 


0.04  - 

I 

0.03  h 

I 

0.02  L 

I 

0.01  ^ 
0.00  p 

-0.01  I 

I 

-0.02  'p 

-0  03  ^ 

1 

-0.04  ' — 
0.0 


AEDC  Wind  Tunnel  Data 
Leading  Edge  Angle 
■  =45  deg.  -TTCP  standard 
•  “180  deg.  -blunt 


■  ■ 


_  ' _ . _ _ _ j__^ _ _ _ 

05  1.0  15  20  2.5  3.0  35 

MACH  NUMBER 


Roll  Moment  CoefRcient  vs.  Mach  Number  for  Configurations 
With  Different  Leading  Edge  B luntness . 


23 


ROLL  MOMENT  COEFFICIENT 
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